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ABSTRACT 

The CMB has distinct peaks in both its temperature angular power spectrum (TT) 
and temperature-polarization cross-power spectrum (TE). From the WMAP data we 
find the first peak in the temperature spectrum at i = 220.1 it 0.8 with an amplitude 
of 74.7 ± 0.5 fiK; the first trough at £ = 411.7 ± 3.5 with an amplitude of 41.0 ± 0.5 ^K; 
and the second peak at £ = 546 it 10 with an amplitude of 48.8 it 0.9 ;uK. The TE 
spectrum has an antipeak at ^ = 137 it 9 with a cross-power of —35 it 9 /uK^, and a peak 
at £ = 329 lb 19 with cross-power 105 it 18^K^. All uncertainties are la and include 
calibration and beam errors. 

An intuition for how the data determine the cosmological parameters may be gained 
by limiting one's attention to a subset of parameters and their effects on the peak 
characteristics. We interpret the peaks in the context of a flat adiabatic ACDM model 
with the goal of showing how the cosmic baryon density, matter density, Q^nh'^, 

scalar index, n^, and age of the universe are encoded in their positions and amplitudes. 
To this end, we introduce a new scaling relation for the TE antipeak-to-peak amplitude 
ratio and recompute known related scaling relations for the TT spectrum in light of 



^ WMAP is the result of a partnership between Princeton University and the NASA Goddard Space Flight Center. 
Scientific guidance is provided by the WMAP Science Team. 

^Dept. of Physics, Princeton University, Jadwin Hall, Princeton, NJ 08544 

^Code 685, Goddard Space Flight Center, Greenbeh, MD 20771 

^Dept. of Physics and Astronomy, University of British Columbia, Vancouver, BC Canada V6T IZl 
^Depts. of Astrophysics and Physics, EFI and CfCP, University of Chicago, Chicago, IL 60637 
^Dept. of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 
^Dept. of Physics, Brown University, Providence, RI 02912 
*UCLA Astronomy, PO Box 951562, Los Angeles, CA 90095-1562 
^National Research Council (NRC) Fellow 



- 2 - 



the WMAP data. From the scaling relations, we show that WMAP^s tight bound on 
Jlfe/i^ is intimately linked to its robust detection of the first and second peaks of the TT 
spectrum. 

Subject headings: cosmic microwave background, cosmology: observations 

1. Introduction 

WMAP has mapped the cosmic microwave background (CMB) temperature anisotropy over 
the full sky with unprecedented accuracy (Bennett et al. 2003b). The temperature angular power 
spectrum (TT, Hinshaw et al. 2003) and the temperature-polarization cross-power spectrum (TE, 
Kogut et al. 2003) derived from those maps have a number of characteristic features. It is these 
features, and our ability to predict them, that make the anisotropy such a powerful tool for cosmol- 
ogy. Computer programs like CMBFAST (Seljak &: Zaldarriaga 1996) efficiently compute the TT 
and TE power spectra for a wide variety of cosmological parameters. The model spectra are then 
compared to the data to deduce the best fit parameters (Spergel et al. 2003). The distinctiveness 
and accuracy of the measured spectrum determines the degree to which the parameters may be 
distinguished. 

While the parameters of cosmological models are ultimately found by maximizing the likelihood 

of the data given a model, and the validity of a model is assessed by the goodness of fit (Spergel 
et al. 2003), this analysis by itself provides no model independent assessment of the features of the 
angular power spectrum. Both intuition and calculational simplicity fall by the wayside. In this 
paper we focus on the peaks of the TT and TE spectra. There are four reasons to consider just 
these particular characteristics. The first is that through an examination of the peaks, one can 
gain an intuition for how the cosmological parameters are encoded in the TT and TE spectra. The 
second is that by determining how the peaks depend on cosmological parameters, one may quickly 
assess how potential systematic errors in the angular power spectra afi^ect some of the cosmological 
parameters. The third is that alternative models, or additions to the best fit model, may be easily 
compared to the data simply by checking to see if the alternative reproduces the peak positions 
and amplitudes. Finally, the peaks serve as a simple check for the sophisticated parameter fitting 
described in Spergel et al. (2003) and Verde et al. (2003). 

The paper is organized as follows. In §2 we give a brief overview of the CMB. In §3 we give 
WMAP^s best fit values for the positions in ^-space and amplitudes of the peaks. We then consider, 
in §4, the cosmological information that comes from the positions and amplitudes of the TT and TE 
peaks. We interpret the peaks in terms of the best fit model to the full WMAP data set (Spergel 
et al. 2003) and make direct connections to the decoupling epoch. We conclude in §5. 
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2. Overview 

The CMB TT spectrum may be divided into three regions depending on the characteristic 
angular scale of its features. In each region, a different physical process dominates. The regions 
correspond to: (a) Angular scales larger than the horizon size at decoupling as observed today. 
These correspond to 9 > 2° or, cquivalcntly i < idee ~ 90. The low i portion is termed the 
Sachs- Wolfe plateau (Sachs &: Wolfe 1967). In this region one observes the relatively unprocessed 
primordial fluctuation spectrum because patches of sky with larger separations could not have been 
in causal contact at decoupling, (b) The acoustic peak region, 0?2 < ^ < 2° or 90 < £ < 900, which 
is described by the physics of a 3000 K plasma of number density ~ 300 cm~^ responding to 
fluctuations in the gravitational potential produced by the dark matter, (c) The Silk damping 
tail (Silk 1968), 6 < 0?2 or i > 900, which is produced by diffusion of the photons from the 
potential fluctuations and the washing out of the net observed fluctuations by the relatively large 
number of hot and cold regions along the line of sight. The basic framework in which to interpret 
the temperature anisotropy has been known for over thirty years (Peebles & Yu 1970; Sunyaev &: 
Zel'dovich 1970). 

It has long been recognized that the positions and amplitudes of the peaks in the 90 < £ < 900 
region could be used to constrain the cosmological parameters in a particular model (Doroshkevich 
et al. 1978; Kamionkowski et al. 1994; Jungman et al. 1996b, a). Previous studies (Scott et al. 1995; 
Hancock & Rocha 1997; Knox k Page 2000; Weinberg 2000; Cornish 2001; Podariu et al. 2001; 
Miller et al. 2002; Odman et al. 2002; Douspis & Ferreira 2002; de Bernardis et al. 2002; Durrer 
et al. 2003) have steadily increased our understanding of the peaks and their significance. 

A summary of pre- WMAP determinations of the position and amplitude of the first peak is 
given in Table 1. 

3. Determination of Peak Characteristics 

We determine the peaks and troughs with fits of Gaussian and parabolic functions to the data. 
Through this process, we compress a large data set to eight numbers. Such a compression is similar 
to specifying the cosmological parameters; though, it is considerably easier to compute and directed 
more toward the intrinsic characteristics of the data as opposed to a cosmological model. 

For the TT spectrum, AT/ = £{£ + l)Q/(27r), we model the spectrum for 100 < ^ < 700 
as composed of a Gaussian peak, a parabolic trough, and a second parabolic peak independent of 
any particular cosmological model. Using a parabola to fit the first peak results in a higher 
and systematically underestimates the amplitude. Parabolas are adequate for the other features. 
The independent parameters are the amplitude and position of each peak and trough, the width 
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of the first peak, and a continuity parameter. The latera recta of the trough and second peak 
are derived by constraining the model spectrum to be continuous. The junction between the first 
peak and the first trough (first trough and second peak) is fixed to be where the first peak (first 
trough) equals the continuity parameter. In other words, the continuity parameter is the value of 
the angular power spectrum where the curves meet. There are a total of eight parameters. We 
denote the positions of the first peak, first trough, and second peak as , iJJ, £2'^ respectively. 
The amplitudes of the first and second peak are AT^j^j, and Ar|yy. 

The peak parameters are found in three ways: (a) with a direct non-linear fit based on the 
Levenberg-Marquardt algorithm (Press et al. 1992); (b) with a Markov Chain Monte Carlo (MCMC) 
(Christensen et al. 2001) of the peak parameters (not the cosmological parameters); and (c) with 
peak-by-peak fits with a Gaussian fitting function. We have also separately fit the first peak with 
a Gaussian with additional parameters for kurtosis and skewness. All methods yield consistent 
results. The values we quote are from the MCMC. All fitting is done with the unbinned angular 
spectra using the full covariance matrix S^^/ (Verde et al. 2003) . The year one WMAP calibration 
uncertainty is 0.5% which is added in quadruture to the fitted peak amplitude uncertainty. The fit 
to the TT peak models yields a minimum x^/^ = 702/593 = 1.18. While not as good a fit to the 
data as a CMBFAST-derived spectrum, the determinations of the peaks are the same for all three 
methods and are consistent with those of the best fit CMBFAST model. 

For the TE spectrum the data between 40 < / < 450 are modeled as a piecewise-continous 
composite of a parabolic antipeak and parabolic peak. (We do not fit the reionization region, 
£ < 20). Each parabola has a latus rectum, height, and position; however, the latus rectum of 
the antipeak is constrained by the requirement that the antipeak and peak are continuously joined 
at the zero-crossing. Thus there are just five free parameters. We impose a prior on the latus 
rectum of the peak to be < 150. The TE Fisher matrix depends on Cf^, Cj^, and Cf^. The TT 
spectrum is fixed to be a best-fit model, and the small EE contribution is neglected. The best fit 
model spectrum has = 504/466 = 1.08. 

The parameters are summarized in Table 2 along with the peak positions determined directly 
from the MCMC method. These two completely distinct methods yield consistent results. Figure 1 
shows the binned TT and TE angular spectrum with the best fit peaks model and the la and 2a 
contours for the peak and trough positions. The amplitude of the TT trough and second peak are 
separated by > 5cr, leaving no doubt of the existence of the second peak. 

Also shown in Table 2 are the peak determinations from the full analysis using just the WMAP 
data (Spergel et al. 2003). The agreement between the two is generally good, most of the values 
are within la, all are within l.Qa. This level of discrepancy is expected because the peak fitting 
and the full analysis weight the data in very different ways. 



The latus rectum of a conic section is the length of the chord parallel to the directix that passes through the 
focus. For the parabolic curve y = yo + {x — xof /w the value of the latus rectum is |w|. 
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4. Interpretation of Peaks and Troughs 

We now interpret the peak characteristics in terms of a flat adiabatic ACDM cosmological 
model. We focus on the baryon density, ui, = flbhP, the matter density, ujm = ^mh"^ {^m = ojjj + uic, 
where Uc is the cold dark matter component), and the the slope of the primordial power spectrum. 
Us- The comparison of thc^ j^caks and troughs to the model is made at two levels. At the more 
general level, we use CMBFAST to derive simple scaling relations and show how the cosmological 
parameters are encoded in the peak characteristics. The aim is to build an intuition for the 
connection between the raw data and the deduced parameters. This approach has been championed 
by Hu et al. (2001); this paper forms the basis for our presentation. The scaling relations hold for 
a wide range of parameters and may be derived without considering WMAP data. However, we 
optimize them for the best fit WMAP parameters (Table 3), so they arc as accurate as possible. 
At a more detailed level, we show that by using the peak values and uncertainties with the scaling 
relations, one obtains constraints and uncertainties on a;;,, uim, ns, and the age of the universe that 
are consistent with the full analysis (Spergel et al. 2003). 

Seljak (1994) and Hu & Sugiyama (1995) showed that the physics of the acoustic peaks may 
be understood in terms cob, ^m-, f^s-, t (the optical depth to when the universe was reionized), and 
6a (the angular scale of the sound horizon at decoupling) . Our approach follows that of Hu et al. 
(2001) in that we consider just the positions of the peaks and the ratios of the peak amplitudes. 
For i > 40, the peak ratios are insensitive to the intrinsic amplitude of the CMB spectrum and to 
r. The free electrons from reionization at z 20 (Kogut et al. 2003) scatter CMB photons, thereby 
reducing the CMB fluctuations by nearly a constant factor for multipoles i > 40. (The reionization 
also increases the TT spectrum at ^ < 20 due to Doppler shifts of the scatterers.) We therefore 
consider just oJb, ujm, rig, and 6a, in the context of a flat ACDM model. 



4.1. The Position of the First Peak 

The acoustic peaks arise from adiabatic compression of the photon-baryon fluid as it falls 
into preexisting wells in the gravitational potential. These potential wells are initially the result 
of fluctuations in some primordial field (e.g., the inflaton field in inflationary cosmology). The 
wells are enhanced by the dark matter (ujc), which is able to cluster following matter-radiation 
equality (z^q, Table 3) because the dark matter, for our purposes, does not scatter off the photons 
or baryons. 

The first peak corresponds to the scale of the mode that has compressed once in the age of the 
universe at the time that the photons decoupled from the electrons at z^ec, some 379^7 kyr after 
the Big Bang (Spergel et al. 2003). The characteristic angular scale of the peaks is set by 

dAiZdec)' 

where is the comoving size of the sound horizon at decoupling and dA is the comoving angular 
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size distance to the decoupling surface. 

The size of the sound horizon may be computed from the properties of a photon-baryon fluid 
in an expanding universe. At decoupling it is given by (Hu &; Sugiyama 1995, B6): 



r,(.,ec) = 3997,/^ In + ^^ec + Vg.ec + ^ 

^ ' y UJmCOb 1 + V^ 

in Mpc, where R{z) = 3pb{z)/4:pj{z) with the baryon density and the photon density. The 
redshift at matter-radiation equality is 

, ^ 5464 (u;^/0.135) 

(rcMB/2.725)4(l + p,/p^)- ^> 

The ratio of neutrinos to photons is Pu/p-y = 0.6851 (Hannestad & Madsen 1995). Note that ^^(zdec) 
depends only on the physical densities, ujm and ujf,, and not on h. (We assume Tcmb = 2.725 K 
(Mather et al. 1999) and the number of neutrino species, N^, = 3.) Also, Vg is independent of the 
curvature and cosmological constant densities, ilk and Oa, because these are both late acting: at 
high z the universe evolves as though it is geometrically fiat. 

The comoving angular size distance^-*^ to the decoupling surface, d^, for a flat geometry is 

= r . (4) 

where Qr is the current radiation density and = 1 — — Qj. (Peebles 1993). Neither nor 
cIa depend on Ug. The picture one has is that of measuring the angular size 9a of a physical meter 
stick of size near the edge of the observable universe with geodesies of the intervening geometry 
as encoded in (Ia- 

The position in ^-space of the first peak is intimately linked to 6a- In fact, in the full analysis 
(Spergel et al. 2003), 6a is one of the best determined fit parameters because of its association with 
the first peak. To make the connection, one defines a characteristic acoustic index, £a = tt/^a- 
However, Ia is not ij'^ because the potential wells that drive the compression, leading to the 
large variance at ij"^, are dynamic: they respond to the matter and radiation that falls into them. 
Additionally, there is no precise physical relation between an angular scale and an associated i. 
One calibrates the relation with cosmological models, such as those from CMBFAST, to find a 
phase factor (pi that relates the two (Hu et al. 2001). The general relation for all peaks and troughs 
is 

C = ^A(m-0^), (5) 



^"'We follow the naming convention of Peebles (1993, 13.29) for d^- Common alternative names include "proper 
motion distance" (e.g., Weinberg 1972, 14.4.21), "angular diameter distance" (e.g., Efstathiou & Bond 1999), and 
"transverse comoving distance" (e.g., Hogg 1999). 
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where m labels the peak number (m = 1 for the first peak, m = 1.5 for the first trough, etc.). For 

a particular curvature density, Q^, the <f>m depend weakly on cuj, and cOm- For example, near the 
WMAP values, changing Ug, ujf,, and by 10% changes 0i by 4.6%, 0.5%, and 1.1%, respectively. 
Because the (prn differ by only « 25% between peaks, £a may be thought of as the characteristic 
scale of the peaks in ^-space. Effectively, the full analysis solves for the phase factors simultaneously 
with the other parameters. 

The considerations above are quite general and involve a good fraction of the quantities that 
are of interest to cosmologists. We now show how the WMAP peak positions are related to, and 
in some cases determine, these quantities. 

The position of the first peak is an essential ingredient of the mixture of properties that 
make the CMB such a powerful probe of the geometry of the universe (Kamionkowski et al. 1994; 
Jungman et al. 1996b). However, the position alone does not determine the geometry. With the full 
TT spectrum, the quantities that WMAP measures particularly well, independent of f^^, are iOm 
and LOh because both these quantities affect the spectrum at early times, before geometric effects 
shift the spectrum. Even if ujm and u>i, are known, il^. and Q,m (through h) may be traded off 
one another to give the same 6a- This is called the "geometric degeneracy" (Bond et al. 1997; 
Zaldarriaga et al. 1997). The degeneracy is broken with either prior knowledge of h or $7^ (e.g., 
Freedman et al. 2001; Bahcall et al. 1999). For WMAP if we impose a prior constraint of /i > 0.5 
then = 0.03 ±0.03; without the prior, = -0.050^°;o39 (Spergel et al. 2003). In the following, 
we take the geometry to be flat and therefore consider the peak positions as a function of Uh, Ug, 
and ujm only. 

With known ui, and the acoustic horizon size is particularly well determined. From 
equation 2 we find = 143 ± 4 Mpc. When Ug is included, we may in addition determine the 
phases. From the fits in Doran & Lilley (2002) we obtain = 0.265 it 0.006. The uncertainty is 
derived from the uncertainties of Hg, ujb, <^rm and the quoted accuracy of the fitting function. This 
combined with the measured position of the first peak at (.f^ = 220.1 it 0.8 gives us an acoustic 
horizon scale of £a = 300 it 3. And, from Ia we find that 9a = 0.601 it 0.005°. In other words, 
we know the angular and physical sizes of structures on the decouping surface very well. From 
equation 1, we solve for the angular size distance and find that dA = 13.7 ± 0.4 Gpc. If we had 
naively tried to compute dA from h and and their uncertainties directly from equation 4, the 
resulting uncertainty would have been larger. This, though, is not the correct procedure because 
the h and flm deduced from the CMB are correlated as shown in Figure 2. 

The interplay between measuring the peak position, and thus dA, and measuring cOm is at 
the root of WMAP''s ability to determine and h separately (Bond et al. 1994; Efstathiou & 
Bond 1999). In particular, Percival et al. (2002) show that 6a, which depends primarily on the 
first peak position, is the same along a line of constant O^/i^ *^. In other words, this is how 
and h scale to keep 6a constant. In Figure 2 we show the constraint from the measured first peak 
position. For each coordinate pair in the plane, we compute ij'^ from equations 1, 2, 4, and 5 
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and then assign the pair a likehhood derived from the measured distribution of £{ = 220.1 ± 0.8. 
(We set LOh = 0.023 for the calculation as discussed below.) The figure also shows the constraints 
from LOrn and the likelihood contours from the full WMAP likelihood analysis. It is evident that 
the departure from a fixed O^/i^ dependence in 9a enables separating and h. The separation 
is by no means complete and constitutes one of the largest parameter degeneracy for WMAP; 
it precludes determining = 1 — 0^ to high precision with CMB data alone. The independent 
determination of and h will improve with time as the uncertainty on Um decreases. Fortunately, 
galaxy surveys are sensitive to ilmh which breaks the degeneracy when the data sets are combined 
(Spergel et al. 2003). 

It is fortuitous that in a flat geometry the position of the first peak is directly correlated with 
the age of the universe (Knox et al. 2001; Hu ct al. 2001). This is seen most easily in the ~ h 
plane in Figure 2. With w = —1, the age depends only on the expansion time scale h'^ and 
and is given by 



(Carroll et al. 1992). We show the isochrons for 12 Gyr through 16 Gyr. One observes that they 
overlap considerably with the position constraint and the degeneracy lines from equation 4. For 
reasonable values of 0,m, the age is seen to be in the neighborhood of 13.6 Gyr. The full analysis 
(Spergel et al. 2003) gives to = 13.6 ±0.2 Gyr for w = —1. Changing uob changes the position 
constraint through changes in and 01. If Ub is increased from 0.024 to 0.025 (Icr), the position 
constraint line shifts by < a/2, indicating that the plot is not particularly sensitive to changes in 

To summarize this section, we have shown how the position of the first peak is related to the 
cosmological parameters. Traditionally, one has used equation 5, the weak dependence of (j) on ujb 
and ujrm ?^s ~ 1 and the measured position of the first peak to deduce flatness. Instead, we assumed 
flatness and use WMAP's values of Ug, uJb and a;„j, to deduce a precise relation for the acoustic 
scale, iA- We then showed how one of the largest parameter degeneracy in the CMB data could be 
understood in terms trading off h and in 6 a- Finally, we showed that WMAP^ measure of the 
age of the universe is largely a function of WMAP''s identification of the position of the first peak. 



The scaling of the amplitude of the first peak with ujb and ojm has a straightforward interpre- 
tation. However, unlike the position, the amplitude itself has a complicated dependence on the 
cosmological parameters (Efstathiou & Bond 1999; Hu et al. 2001). We defer discussing the Ub 
scaling until after §4.3 and consider here just the ujm scaling. 

Increasing decreases the first peak height. When the universe is radiation dominated at 
z > Zeq and the photon-baryon fluid compresses in a gravitational potential well, the depth of the 




(6) 



4.2. The Amplitude of the First Peak 
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well is decreased due to the additional mass loading of the fluid. As a result, the compressed fluid 

can expand more easily because it sees a shallower potential. The process is termed "radiation 
driving" (Hu k, Sugiyama 1995; Hu &: Dodclson 2002) and enhances the oscillations. (In addition, 
a shallower potential well means that the photons are redshifted less as they climb out enhancing 
the effect further.) The effect remains important through recombination. Increasing u)^, while 
holding all else but f^A constant, moves the epoch of matter-radiation equality to higher redshifts 
without significantly affecting Zc^ec- This gives the dark matter more time to develop the potential 
wells into which the photon-baryon fluid falls. With better defined potential wells, the compressed 
photon-baryon fluid has less of an effect on them, thereby reducing the effects of radiation driving. 
As a consequence, increasing uj^n decreases the amplitude of the fluctuations for all the acoustic 
peaks. 



4.3. Additional TT Peaks and Troughs 

There is a pronounced second peak in the WMAP TT spectrum at a height > 5(7 above the 
first trough. The peak arises from the rarefaction phase of an acoustic wave. In broad terms (Hu & 
Dodelson 2002), in the same conformal time that it takes the plasma to compress over the acoustic 
horizon, an acoustic wave with half the wavelength (twice the i) of the first peak can compress and 
rarify. Likewise, a compressional third peak is expected at the second harmonic of the first peak 
and a rarefaction fourth peak is expected at the third harmonic of the first peak. This set of peaks, 
the first two of which WMAP clearly sees, shows that the photon-baryon fluid underwent acoustic 
oscillations. 

The position of the first trough, and the second and subsequent peaks is predicted from the 
acoustic scale and phase shifts. From equation (5) and the values in Table 3, one predicts l^T, = 
409 ± 4, e'^'^ = 533 ± 5, and ij^ = 809 ± 7. The first two of these values agree with those found 
from the peak fits to within the measurement uncertainty as shown in Table 2. 

The amplitude of the second peak depends on many of the same parameters as the amplitude 

of the first. Its most distinctive feature, though, is that increasing tob decreases its height. The 
reason is that as one increases the baryon density one increases the inertia in the photon-baryon 
fluid. When compared to the lower cob case, the compressions are deeper and the rarefactions not as 
pronounced. Thus, the compressional peaks are hotter (first and third) and the rarefaction peaks 
(second and fourth) are cooler. 

Following Hu et al. (2001) we characterize the amplitude of second peak as a ratio to the 
amplitude of the first peak, HJ^. The ratio is insensitive to the reionization optical depth and 
to the overall amplitude of the spectrum since they simply scale the amplitudes of both peaks. 
It depends on just Ub, tOm, and Ug- The dependence on Um is relatively weak, since to a first 
approximation it too just scales the peak amplitudes as discussed above. The dependence on Us 
simply comes from the overall slope of the CMB angular spectrum. Increasing Us increases the 
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height of the second peak relative to the first. The following function is derived from fitting to 
a grid of CMBFAST spectra. It gives the ratio of the peak amplitudes to 2% accuracy for 50% 
variations in the parameters: 

^TT^J^A^y ^ 0.0264a;,-0-^62(2.42)n.-i 



^g-0.4761n(25.5u)6+l-84u;^)2 



the parameter dependences are 



= 0.88An, - 0.67^ + 0.039^. (8) 

Thus a 1% increase in both U}, and with Ug fixed reduces the height of the second peak relative 
to the first by 0.63%. 

For the WMAP data, i^f"^ = 0.426 ± 0.015. In Figure 3 we show the constraints from //J"^ for 
= 0.14. The contour lines correspond to the full analysis (Spergel et al. 2003). One sees that 
the two analyses are consistent, though the error from just considering .ff is larger. However, it is 
clear that the uncertainty in for U}, at fixed and u;^ leads to nearly the same uncertainty as 
deduced from the full analysis. Thus, we may interpret WMAP^s ability to determine uj}, as coming 
primarily from the precise measurement of the ratio of the first two peaks. The first two terms on 
the right side of the above equation, as shown in Figure 3, also quantify the — Ug degeneracy 
(Spergel et al. 2003). 

The height of the third peak increases as increases, as discussed above, and increasing 
decreases its height. Thus the ratio to the first peak is not as distinctive as for the second peak in 
terms of ujb and u)m but the long I baseline makes the ratio more sensitive to Ug- Hu et al. (2001) 
give: 

HV ^ {^)^ (0) 

2.17a;^59 3.6"--^ 



[1 + (cc;(,/0.044)2][l + 1.63(1 - u;6/0.071)a;^] 

with parameter dependencies 

= 1.28An, - 0.39^ + 0.46^. (10) 

With Ug fixed, increasing ujb and Um by 1% increases the height of the third peak relative to the 
first by only 0.07%. Measuring the third peak helps mostly in measuring Ug, thereby breaking the 
Ub — Ug degeneracy shown in the left side of Figure 3. For WMAP parameters, equation 10 is 
accurate to 1%. 
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WMAP does not yet clearly measure the third peak. For its height, we use the value of the 

Wang et al. (2002) compilation because it is the most recent and includes calibration error. At 
I = 801, (Ar7'^)2 = 2322 ± UOfiK^. Though not quite at the peak, the value is sufficient. With 
WMAP's first peak, Hj^ = 0.42 ± 0.08. 



4.4. The Temperature-Polarization peaks 

The E-mode polarization in the CMB arises from Thompson scattering of a quadrupolar ra- 
diation pattern at the surface of last scattering. At decoupling, the quadrupolar pattern is pro- 
duced by velocity gradients in the plasma and is correlated with the temperature anisotropy (Bond 
Sz Efstathiou 1984; Coulson ct al. 1994). In rough terms, because it is velocity dependent, the 
temperature-polarization correlation traces the derivative of the temperature spectrum. For a 
given TT spectrum, the (non-reionized) TE and EE spectra are predicted. Thus, their observation 
confirms that the cosmological model is correct but in so far as they are predicted for i > 40, they 
contain little additional information about the parameters. By contrast, the TE polarization signal 
at £ < 20 (Kogut et al. 2003), produced by reionized electrons scattering the local Zr ~ 20 CMB 
quadrupole, breaks a number of parameter degeneracies and considerably enhances the ability to 
extract the cosmic parameters. 

From parameterizing the output of CMBFAST, we find that the ratio of the first TE antipeak 
to the second TE peak is given by 

Hr = - (^)' = 0.706a;r'(0.518)--i 

^g0.1951n(33.6a;(,+5.94a;^)2_ ^-|^-^^ 

The parameter dependencies are 

— ^ = -0.66An3 + 0.095 ^+0.45 -. (12) 

The sign preceding Aug is opposite to that in the case of the TT peaks since we are considering 
the ratio of the first antipeak (low £) to the second peak (high i), whereas in the case of the TT 
peaks we were considering the second peak (high i) to the first peak (low i). Here, increasing 
Um increases the contrast between the TE peak and antipeak. This occurs because increasing cOm 
leads to deeper potential wells at decoupling. In turn, this produces higher velocities and thus an 
enhanced polarization signal. From the data we find Hj^ = 0.33 ± 0.10. 



4.5. Combined Peak Constraints 



We now use the scaling relations just presented in a flat model with the scalar index limited 
to Ug = 0.99. This simplifled minimal model leads to a greater intuition for assessing the peaks' 
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role in determining the cosmological parameters. The constraints from H2 and are shown in 
Figure 3. We see that determines cu^ for most values of the dark matter and that wj, < 0.029 

(2cj) regardless of the amount of dark matter. The width of the swath is about twice that of the 
formal error, again indicating that more in the spectrum than just i^^"^ constrains coi,. 

The constraints from H2 and Hi^ arc almost orthogonal to the H2 constraint but their 
uncertainty bands are wide. This shows us that the WMAP value for Uc (or a;^) is not being 
driven specifically by the peak ratios. Rather, the constraint comes from the amplitude of the 
whole spectrum. The degeneracy with r is broken by the TE detection of reionization and the 
degeneracy with the overall amplitude is broken because Um affects primarily the acoustic peak 
region more than the low I part of the spectrum. 

5. Discussion 

We have focused on the dominant acoustic features in the WMAP spectrum, but there are 
other features as well. Notable by its absence is a trough in the TT spectrum at ^ ~ 10 due to the 
increase in fluctuation power near 1 = 1 from the integrated Sachs- Wolfe effect in a A-dominated 
universe. This feature is largely obscured by cosmic variance. For WMAP though, there is very 
little power at low i and the correlations with d > 60° are unusually small (Bennett et al. 2003a). 
The lack of even a hint of an upturn at low i, in the context of ACDM models, is surprising. 
(However, COBE also found no evidence for such an upturn). This departure from the model 
constitutes a small (though possibly important) fraction of the total fluctuation power and does 
not cast doubt on the interpretation of the £ > 40 spectrum. 

At £ ~ 40 and i ~ 210, for example, there are excursions from a smooth spectrum that are 
somewhat larger than expected statistically. While intriguing, they may result from a combination 
of cosmic variance, subdominant astrophysical processes, and small effects from approximations 
made for the year one data analysis (Hinshaw et al. 2003). At present, we consider them interesting 
but do not attach cosmological significance to them (e.g., Peiris et al. 2003; Bose & Grishchuk 2002). 
More integration time and more detailed analyses are needed to understand how they should be 
interpreted. 

In summary, the characteristics of the peaks of the WMAP angular power spectra may be ro- 
bustly extracted from the data. The second TT peak is seen with high accuracy; the TE antipeak 
and peak have been observed for the first time. In the context of a flat adiabatic ACDM model, 
which fits the WMAP data well, we report the characteristics of the decoupling surface. By consid- 
ering a reduced parameter space, we show how the position of the first peak leads to WMAP's tight 
determination of the age of the universe and how ujj, is determined primarily from the amplitude 
ratio of the first to second TT peak. The determination of LOm comes from considering the full data 
set and is not attributable to any particular peak ratio. 
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Table 1. Previous Measurements of the First Peak 



Experiment 



ATia'T [/iK] Cal. error 



£ range 



Reference 



TOCO 




87tl 


8% 


60<e< 410 


Miller ct al. (1999) 


BOOMERANG-NA 


2071^6 




8% 


50 < £ < 410 


Mauskopf at al. (2000) 


MAXIMA 


236 ± 14 




4% 


75 < ^ < 385 


Lee et al. (2001) 


DASI 


19911^ 


71.9 ±4.4 


4% 


110 < ^ < 380 


Halverson ct al. (2002) 


VSA 


2361^0 


73.9 ±6.8 


3.5% 


140 < £ < 420 


Grainge et al. (2003) 


BOOMERANG 


219 ±5 


73.8t?;^ 


10% 


100 < ^ < 410 


Ruhl et al. (2003) 


ARCIIEOPS 


22()t; 




7% 


18 < / < :-!:-!() 


B<-U()it. ct. al. cum) 



WMAP 



220.1 ±0.8 



74.7 ±0.5 



0.5% 



100 < £ < 350 



This paper 



Note. — All values come from fitting a Gaussian shape to just the first peak of the data set specified, 
and include calibration error. Each data set is considered on its own, without the COBE/DMR data, and 
so a direct comparison between experiments may be made. The best fit position depends somewhat on the 
fitting function so the values from different analyses yield different results (e.g., Knox & Page 2000; Durrer 
et al. 2003; Odman et al. 2002; Grainge et al. 2003). The TOCO, VSA, and BOOMERANG-NA experiments 
were calibrated with Jupiter. The TOCO and VSA experiments are most affected because of they operate at 
30-150 GHz and 35 GHz respectively. With the new calibration of Jupiter (Page et al. 2003), the peak values 
above will be reduced « 5%. The weighted peak amplitude is 71.7 ± 2.4 iJ.K and the weighted peak position 
is 218.8 ± 3.5 in good agreement with WMAP. In a separate analysis based on different assumptions, Bond 
reports iV = 222 ± 3 (private communication) . This was also the value preferred by a concordance model 
(Wang et al. 2000) that predated all the experiments of the new millennium. Note that WMAP's values for 
the position and amplitude are both more than four times more precise than all the listed measurements 
combined. 
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Table 2. WMAP Peak and Trough Amplitudes and Positions 



Quantity 


Symbol 


I 




ATI [/iK2] 


FULL^ 


FULL AT| 




First TT Peak 


er 


220.1 ±0.8 


74.7 ±0.5 


5583 ± 73 


219.8 ±0.9 


5617 ± 


72 


First TT Trough 


<^1.5 


411.7±3.5 


41.0 ±0.5 


1679 ± 43 


410.0 ± 1.6 


1647 ± 


33 


Second TT Peak 


pTT 


546 ± 10 


48.8 ±0.9 


2381 ±83 


535 ±2 


2523 ± 


49 


First TE Antipeak 




137 ±9 




-35 ±9 


151.2 ± 1.4 


-45 ± 


2 


Second TE Peak 


pTE 


329 ± 19 




105 ± 18 


308.5 ± 1.3 


117± 


2 



Note. — The values and uncertainties are the maximum and width of the a posteriori distribution of the 
likelihood assuming a uniform prior. The uncertainties include calibration uncertainty and cosmic variance. 
The FULL values are derived from the full CMBFAST-based likeUhood analysis using just the WMAP data 
(Spergel et al. 2003). The FULL method yields consistent results. Recall that the FULL chains are sensitive 
to the combined TT and TE spectra and not just the individual peak regions. Numerical errors in CMBFAST 
will increase the uncertainties, but should not bias the results. 
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Table 3. WMAP Cosmological Parameters for the Peaks Analysis 



Quantity 


Symbol 


Value 




Physical baryon density 




0.024 ±0.001 




Physical mass density 




0.14 ±0.02 


INPUT 


Scalar index 


ria 


0.99 ±0.04 




First TT peak phase shift 


01 


0.265 ±0.006 




First TT trough phase shift 


4>-i-5 


0.133 ±0.007 




Second TT peak phase shift 


02 


0.219 ±0.008 


DERIVED 


Third TT peak phase shift 


03 


0.299 ±0.005 


FROM 


Redshift at dccouphng 


^dec 


1088+2 


INPUT 


Redshift at matter radiation equality 


^eq 


3213I328 




Comoving acoustic horizon size at decoupling (Mpc) 


Ts 


143 ±4 




Acoustic scale 

Comoving angular size distance to decoupling (Gpc) 


Ia 
dA 


300 ±3 
13.7 ±0.4 


DERIVED FROM 
INPUT + PEAKS 



Note. — The cosmological parameters in the top section are derived from just the WMAP data assuming a 
flat ACDM model (Spergel et al. 2003). The quantities in the middle section are derived from the cosmological 
parameters in the top section. The quantities in bottom section are calculated using the middle quantities 
and the measured position of the first peak. The quantity z^ec which wc use corresponds to the location of 
the maximum of the visibility function in CMBFAST. The quantity computed using Hu & Sugiyama (1996) 
corresponds to T{zdec) = 1 and is 1090 ± 2. 



100 200 300 400 500 600 700 
I at Peak and Troughs 



Fig. 1. — The binned WMAP data is shown in blue, the maximum hkehhood peak model from the 
peak fitting functions in red, and the uncertainty contours in black. The top panel shows the TT 
angular power spectrum. The bottom panel shows the TE angular cross-power spectrum. For each 
peak or trough, the contours from the MCMC chains are multiplied by a uniform prior and so they 
are equal to contours of the a posteriori likelihood of the data given the model. The contours are 
drawn at Ax^ = 2.3 and 6.18 corresponding to la and 2a. 
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Fig. 2. — The WMAP data in the - h plane. The thick sohd contours in black are at Ax^ = -2.3, -6.18 
(1(7, 2a) of the marginalized likelihood from the full analysis (Spergel et al. 2003). The filled region is the 
constraint from the position of the first peak, with Wf, = 0.023 fixed. In effect, it shows how rim and h 
must be related to match the observed position of the first peak in a flat geometry, or equivalently to match 
the measured values of 9a- Blue shows the la region and green shows the 2a region. The dotted lines 
are isochrons separated by 1 Gyr. It is clear that the WMAP data pick out 13.6 Gyr for the age of the 
universe in the flat, w = —1 case. The dashed lines show the la limits on The dashed yellow line shows 
^mh^'^ — const. 
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Fig. 3. — Left: Parameter restrictions from H2 iii the uJb — Ug plane. The orange swath is the 
la band corresponding to = 0.426 it 0.015 with uj^ = 0.14. The swath is broadened if one 

includes the uncertainty in Um- The light orange swath is 2a. The solid line in the middle of the 
swath is for AHJ^^ = AuJm = 0. The green contours are from the full analysis of just the WMAP 
data and are thus more restrictive. Right: The constraints in the coh — ujc plane from the peak 
ratios in a flat geometry with = 0.99. The dark shaded regions in each swath are the la allowed 
range; the light shaded regions show the 2a range. Orange is for i^J"^ = 0.426 it 0.015, blue is 
for H'^'^ = 0.42 ± 0.08, and red is for = 0.33 ± 0.10. The uncertainty band for HJ^ is not 

shown as it is broader than the ^^J"^ swath. The heavier central lines correspond to AffJ^ = 0, 
AH^^ = 0, and AHJ^^ = 0, each with Arig = 0. As the mission progresses, all uncertainties will 
shrink. 



